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We prove that a homogeneous and isotropic universe containing a scalar field with a power-law 
potential, V (0) = AcfP, with 0 < n < 1 and A > 0 always develops a finite-time singularity at which 
the Hubble rate and its first derivative are finite, but its second derivative diverges. These are the 
first examples of cosmological models with realistic matter sources that possess weak singularities of 
’sudden’ type. We also show that a large class of models with even weaker singularities exist for non¬ 
integer n > 1. More precisely, iffc<n<fc-|-l where fc is a positive integer then the first divergence 
of the Hubble rate occurs with its {k + 2)th derivative. At early times these models behave like 
standard large-field inflation models but they encounter a singular end-state when inflation ends. 
We term this singular inflation. 

PACS numbers: 98.80.Bp, 98.80.Cq, 98.80.Es, 98.80.Jk 


I. INTRODUCTION 

When can a universe possess a future singularity? [1] 
In the absence of quantum effects, the best known case 
where this can happen, first found by Friedmann [2], is 
a future ’big crunch’ singularity of infinite density. It is 
essentially the time reversal of the big bang singularity: 
the universe expands to a maximum size at some time, 
contracts, and then collapses to the future singularity. 
This occurs in an isotropic universe if the cosmological 
constant, A, is negative, or if the universe is positively 
curved and contains a perfect fluid with bounded pressure 
obeying the strong energy condition. 

A much more drastic type of future singularity is the 
’big rip’. This was first studied in Refs. [3-5]. It oc¬ 
curs when the universe contains a fluid violating the null 
energy condition. The expansion is so rapid that the ex¬ 
pansion scale factor and all its time-derivatives diverge 
in finite time. For a perfect fluid with equation of state 
p = wp, linking its pressure p and density p, a ’big rip’ 
occurs when w < —1. 

By contrast, a far weaker type of finite-time singularity 
is a ’sudden singularity’. Strikingly, it does not require 
cosmological contraction to occur. These singularities 
were first discovered in Ref. [6] as a counterexample to 
the belief (Ref. [7]) that compliance with the strong and 
weak energy conditions and the positive pressure crite¬ 
rion alone suffice to ensure that a spatially closed uni¬ 
verse with topology has a maximal hypersurface and 
recollapses. They were later discussed systematically in 
Refs. [8, 9] (see Refs. [10, 11] for similar singularities 
in brane-world models). They are characterised by the 
scale factor, a{t), and its first derivative d{t) being finite 
as t —>■ ts, but d{t) Qo as t ^ tg. At these singularities 
p is finite but p -boo, and p ip > 0. Generalised 

sudden singularities can also be constructed where the 
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singularity occurs in arbitrarily high derivatives of a{t), 
with all lower derivatives finite [9]. For further studies 
of these classical singularities see Refs. [12-18] and for 
discussions of quantum effects see Refs. [19, 20]. 

Sudden singularities stand in stark contrast to big rip 
singularities as they are much weaker. Firstly, they sat¬ 
isfy all of the classical energy conditions bar the domi¬ 
nant energy condition (and generalised sudden singulari¬ 
ties satisfy them all) [21]. Secondly, while they are scalar 
polynomial curvature singularities, they are ’weak’ in the 
sense of both Tipler [22] and Krolak [23] (this means they 
are not “crushing” singularities because an object ap¬ 
proaching the singularity is not crushed to zero volume). 
Moreover, geodesics are actually extendible through sud¬ 
den singularities [24], unlike for big rip and big crunch 
singularities^, and this behaviour is stable [25]. For other 
types of cosmological singularities discussed in the liter¬ 
ature see Refs. [26-31]. 

Unfortunately, unlike the known examples of a big 
rip singularity, we do not have a simple physically well- 
motivated matter model which produces a sudden (or 
generalised sudden) singularity. Most studies of these 
singularities postulate a form of the scale factor with the 
desired singularity and then use the Einstein equations 
to find energy-condition compliant p and p behaviours 
which source this solution - although a solution of this 
sort is general in the function-counting sense [14]. 

Here, we will construct for the first time a large fam¬ 
ily of finite-time cosmological singularities for a canonical 
scalar field with a simple, power-law potential. We show 
that the formation of these singularities is a general fea¬ 
ture of these simple matter models. They provide the 


^ Note in this paper the term ’singularity’ is used somewhat loosely 
to denote a point where some, in principle, observable quantity 
(such as H or its higher order derivatives) becomes unbounded at 
finite time. It is not necessarily accompanied by geodesic incom¬ 
pleteness, and thereby the Hawking-Penrose singularity theorems 
do not apply in this case. The somewhat informal terminology is 
common in the literature in this subject; a better term for these 
very weak singularities might be singular events. 
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first examples of cosmological models with realistic mat¬ 
ter which evolve towards weak, finite-time singularities. 
They also have the interesting feature that, for appro¬ 
priate initial conditions, they can describe large-field in¬ 
flation in the early universe. The qualitative difference 
to standard inflation models is that when inflation ends 
they evolve to a singular state in finite time. 

II. SCALAR FIELD COSMOLOGIES 

We consider a spatially-flat Friedmann-Lemaitre- 
Robertson-Walker (FLRW) universe, with scale factor 
a{t) and Hubble parameter H{t) = a{t)/a{t). We assume 
the universe contains only a scalar field, (pit), with self¬ 
interaction potential V{(j)). The Einstein and scalar field 


equations of motion are (in units where c = 8ttG = 1) 

3H^ = ^P^ + V{P), ( 1 ) 

H = ( 2 ) 

P + 3HP + V’{(I)) =0, (3) 

where ' = d/dp. We consider the case where the potential 
takes a power-law form: 

V{P) = AP-, (4) 


where A > 0 is a constant. When n is a positive even 
integer it provides the classic example of a large-field in¬ 
flation model with a single potential minimum. If the 
initial conditions are chosen so that the scalar field starts 
high enough up the potential then the universe inflates 
as the p field rolls slowly down the potential. When it 
nears the potential minimum at </> = 0 inflation ends and 
the field oscillates about the minimum. In particular, the 
case where n = 4 is the original model of chaotic inflation 
[32]. When n is a positive odd integer the universe ap¬ 
pears to recollapses under the influence of the scalar field 
(for the n = 1 case see Ref. [33]). We will be interested 
in the case where n > 0 and is not an integer. 

A. Finite-time singularities when 0 < n < 1 

We first examine the case where 0 < n < I. Eq. (3) 
becomes 

P = -3HP- AnP'^-^. (5) 

At t = 0 we choose initial conditions so that po = (/)(0) > 
0 (cases where po < 0 are not cosmologically relevant, 
and generally unphysical), but the value of po = p(0) is 
unconstrained. We assume that the universe is expanding 
initially, so Hq = i^(0) > 0. It is not difficult to see how 
the system evolves in time. If we start with po > 0 
then both terms on the right-hand side of Eq. (5) are 
negative, so in a finite time p becomes negative. Hence, 


in finite time the scalar field starts to decrease, and since 
p continues to decrease, because the second term on the 
right-hand side of Eq. (5) increases as p decreases, it will 
reach <)) = 0 in finite time. When this happens p ^ —oo 
(since n — I < 0) as <()—>• 0. From Eqs. (I)-(2), we see H 
and H are both finite at this point (provided that p <Q 
is finite), but H diverges as 

H = —pp —>■ —oo as p ^ 0. (6) 

Before we demonstrate these claims rigorously we make 
a few points about the nature of the singularity encoun¬ 
tered. Firstly, note it is not a scalar polynomial (sp) 
curvature singularity, as both H and H are finite at this 
point. For a spatially-flat FLRW spacetime the Ricci 
scalar, R, may be written as 

R = 6{2H^ + H), (7) 

which is clearly finite as ^ 0. However, higher scalar 

derivatives of the curvature (like daRd°‘R or Di?) are not 
regular since 

R = Q{AHH + H)^-oo as </>-)> 0. (8) 

These are generalised sudden singularities and examples 
of what Ellis and Schmidt call scalar curvature sin¬ 
gularities [34]. We might call them scalar differential 
singularities. They are the first examples of such weak 
singularities in an FLRW spacetime for a simple mat¬ 
ter model (see Ref. [35] for similar singularities in tilted 
Bianchi spacetimes). Regularity of the curvature ensures 
that they are weak in the sense of both Tipler [22] and 
Krolak [23]. Furthermore, since they are weaker than 
sudden singularities the spacetime is geodesically com¬ 
plete and extendible at the finite-time singularities. 

We can ask which of the classical energy conditions 
these singularities obey. A canonical scalar field is equiv¬ 
alent to a perfect fluid with density p = \(jP + V and 
pressure p = ^p^ — V. This implies that the null energy 
condition {p + p > 0) is always satisfied for any choice of 
V{p). Moreover, if V{p) > 0, as in our case, then both 
the weak energy condition (p > 0 and p + p > 0) and 
the dominant energy condition (p > 0 and p > ]p|) are 
satisfied. The strong energy condition {p + p > 0 and 
P + 3p > 0) is satisfied if 0^ > H; while this initially may 
not hold depending on the initial conditions we choose 
(indeed, if this model is to function as inflation at early 
times it must be violated), it is always satisfied as ^ 0. 
Therefore all classical energy conditions are satisfied near 
the singularity. Again, this is the first example of this 
kind known for a simple matter model. 

Let us now prove that the system does indeed develop 
a finite-time singularity of this type. It suffices to prove 
that the system reaches </> = 0 in finite time with p regu¬ 
lar at this point and ^ < 0. To do so, we prove first that 
the system always reaches (/) = 0 in finite time, and then 
show that it must subsequently reach (^ = 0 in finite time. 
Note that without loss of generality we may assume that 
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00 > 0 , as otherwise we simply omit the first step in the 
argument. 

To prove the first claim, note that in this region of 
parameter space 0 < 0 , so that 0 is strictly decreasing 
and so, while 0 increases in this region, it cannot grow 
faster than linearly in time: 0(t) < 0o + 0oi- In Eq- (5) 
this implies that 

0 < —nA0"“^ < —nA{(j)o + 0o^)"~^- (9) 

Integrating Eq. (9) once gives an upper bound on 0: 

0( i ) < ^ _ (^0 + (10) 

00 00 

This implies that 0 = 0 is reached in finite time. This 
completes the first part of the argument. 

For the second part, we must show that the system will 
reach 0 = 0 in finite time. Before we do so, first note that 
at the point where 0 = 0 we have 0 = —AncfA~^ < 0, so 
there exists a later time, T, at which (j){T) < 0, and 0 
cannot become positive again. We will prove that 0 = 0 
is reached in finite time by showing that 0 is a decreasing 
function of time, by which we mean that 0(t) < 0(T') for 
all t > T. Since 0 is negative this clearly suffices to show 
0 = 0 is reached, as integrating this inequality once gives 
0(i) < 4’{T) + (< —T)0(T), and so 0 = 0 is reached within 
a time t = 0(T)/|0(T)| from t = T. 

To prove that 0 is decreasing, note that without loss 
of generality we may choose T so that 0(T) < 0 as well. 
Now if 0 increases in the future there must exist a time 

> T at which 0(ti) = 0(T) and 0(ti) > 0. However, 
this is impossible since Eq. (2) implies H is a decreasing 
function of time and at this point Eq. (5) requires 

0(ti) = -3i?(ti)0(ti) - Hn0(ti)"“^ 

< -3iI(T)0(T) - Hn0(T)”-i = 0(r) < 0, (11) 

hence no such time exists and so 0 is always decreasing. 

While we have proven that 0 = 0 is always reached 
in a finite time, regardless of the initial conditions, we 
also need to show that 0 is regular at 0 = 0. This is 
easily shown as H{t) is a decreasing function of time and 
positive in this interval, so 


reached in finite time, as from Eq. (5) we have that 
0 < —nA(/>Q~^, which implies that 

0(i) < 00 - nH0p“4, (13) 

so 0 = 0 is reached within time to = 0o/^^0o~^- More¬ 
over, at this point 0 < 0, so 0 becomes strictly negative 
and can never become positive, so 0 always decreases 
from this point. It is not so easy to show analytically 
that 0 = 0 is reached in finite time because the key step 
used in the argument when 0 < n < 1 - that 0 is always 
a decreasing function - is no longer true when n > 1. 
However, it is easy to see that the only alternative is 
that 0 —>■ 0 and 0 ^ 0 as t —>■ oo. This situation is at 
most of measure zero, and it is easy to show via numerical 
simulations of Eqs. (l)-(3) that 0 = 0 is indeed reached 
in finite time for generic initial conditions for any chosen 
non-integer value of n > 1. 

When 1 < n < 2 we can show that 0 = 0 is reached in 
finite time by the following argument. Firstly, without 
loss of generality we may assume that for t > ti > to 
we have 0 > 0, otherwise 0 decreases even faster and so 
0 = 0 must be reached in finite time if it is reached when 
0 > 0. Now 0 > 0 implies that for t > ti we have that 


^ < 


—Ancj)'^ ^ ^ —Ancj)"' ^ 
3ff - SH{ti) 


(14) 


since H{t) is a decreasing function of time. Integrating 
this once gives 


02-n 


An{2 — n)t 

5H{ti) 


(15) 


where C > 0 is a constant. Clearly then, if n < 2 this 
implies that 0 = 0 is reached in finite time. 

When n > 2 one can use the results of Richard [36] to 
derive conditions under which 0 = 0 is reached in finite 
time. Richard studied the existence of zeros of a class of 
second-order, non-linear ordinary differential equations 
which include Eq. (3). Applying his results to Eq. (3) 
we find a zero exists in finite time provided the function 

2 

n + 1 

( 16 ) 


m = -<.(!)” 


{tydf 


JQ 


02(0 = 0) <2A0”(0 = O). (12) 

Since 0 is finite at the turnover we conclude that 0 is 
finite at 0 = 0, and as it is a decreasing function of time 
we must also have (p < 0. 


B. Finite-time singularities when n > 1 

Let us now examine the case where n > 1 and n is 
not an integer. Here, we also expect that 0 = 0 should 
be reached in a finite time from generic initial data. In 
particular, it is possible to show that 0 = 0 is always 


is decreasing and 0 —)> 0 as t oo. Notice that in an 
expanding universe 9{t) > 0, so in our case 0(t) always 
approaches zero from above. In general, 9{t) will be a 
decreasing fnnction provided that a{t) increases slowly 
enough. For instance, for a power-law expansion, a(t) oc 
t™, then 9{t) is a decreasing function which tends to zero 
asymptotically provided that 


m < 


1 

3 


4 

3(n — 1) 


(17) 


That is, provided the expansion rate eventually ap¬ 
proaches that of a kinetically-dominated solution the con¬ 
ditions are satisfied and 0 = 0 is reached in finite time. 
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When (j> = 0 is reached Eq. (12) still implies (()< 0 is fi¬ 
nite. Moreover, Eq. (3) now implies that (p is finite at this 
point, unlike in the previous case. However, higher-order 
derivatives of (p diverge. For instance, differentiating Eq. 

(3) once gives that 

"(P - 9H^<p - ^<p^ - 3HV'{p) + V"{p)<p = 0. (18) 

For 1 < n < 2, every term is finite as ^ ^ 0 except the 
last term which is divergent, so (p ^ oo as (/)—>■ 0. This 
implies that the first divergence in the scale factor occurs 
at fourth order in its derivatives, since 

ii =—ijp'— p 'cp ^ oo as (()—>• 0. (19) 

This implies, for instance, that Di? and higher deriva¬ 
tives of the curvature are divergent on approach to this 
singularity. 

It is not difficult to generalise these conclusions to ar¬ 
bitrarily large non-integer values of n. If fc < n < fc -I- 1, 
where fe is a positive integer, then as (/)—>■ 0 we have 
(jP.k+ 2 ) (—l)^+ic)o, with all lower derivatives of p finite. 

This implies that the first divergence of the Hubble rate 
occurs for the (fc -I- 2)*^ derivative: —>• (—l)^+^oo 

as ^ 0. If n is an integer these singularities never occur 
as V{p) is smooth at ^ = 0. 

III. DISCUSSION 

We have shown that a power-law potential given by Eq. 

(4) admits an arbitrarily large family of ultra-weak gen¬ 
eralised sudden singularities which satisfy all the classical 
energy conditions and are characterised by the divergence 
of a sufficiently high derivative of a curvature invariant. 
Polynomial curvature invariants are always finite upon 
approach to the singularity: the divergence always oc¬ 
curs in a derivative of a curvature invariant due to the 
non-analytic behaviour of the curvature at the singular¬ 
ity. More generally, we expect that any potential V{p) 
which is not smooth at ^ = 0 will admit singularities of 
a similar form. 

It is important to note that these models describe in¬ 
flationary cosmologies in the same way as large-held in¬ 
flation models: we simply choose our initial conditions, 
pQ and pQ, so that the system starts high enough up the 
potential and is potential-dominated. Inhation then oc¬ 
curs while the held violates the strong energy condition, 
but it eventually ends as </>—>■ 0 and the system enters 
the reheating phase. The only difference is when the 
system reaches ^ = 0 deep in the reheating phase (us¬ 
ing the slow-roll approximation this is reached in time 

t « —7-r— —)• The difference is that reheat- 
ing presumably proceeds differently. Since predictions 
for the power spectrum are insensitive to the behaviour 
at reheating, these models will give the same predictions 
for CMB observables as large-held inhation models. In¬ 
hation with fractional potentials has been studied before 


[37], and can even be motivated from string theory [38], 
although their singular behaviour as ^ ^ 0 had not been 
recognised. Note though that the arguments used to de¬ 
rive fractional potentials from string theoretic construc¬ 
tions are only valid for large held values in general. It is 
also worth noting that monomial potentials with n < 2 
given a better ht to current CMB data than large integer 
values do [39]. 

Let us briehy comment on the case when n < 0 and n 
is fractional. These potentials are more relevant as can¬ 
didates for late-time quintessence than inhation. In this 
case, there is little diherence between when n is fractional 
and when it is an integer. One can easily show that if we 
start from initial data where > 0 then, generically, p 
becomes positive in hnite time, and always remains posi¬ 
tive from this point onwards. This is because in this case 
p > 0 for non-zero (() at <(> = 0, so if <() is initially increas¬ 
ing it continues to do so forever. Even if we start from 
initial data where po < 0 then p > 0 and so, by similar 
arguments to section H, p becomes positive in hnite time. 
The scalar held therefore always increases at late times, 
and the evolution is non-singular towards the future. 

To conclude, we have demonstrated the existence of a 
large family of new ultra-weak generalised sudden singu¬ 
larities in a spatially-hat FLRW universe with a scalar 
held possessing a simple power-law potential. They are 
the hrst examples of a spacetime possessing such weak 
singularities for a simple and realistic matter model. 
Moreover, their formation is completely generic in FLRW 
models (indeed, when 0 < n < I they form from any 
isotropic and homogeneous initial data). 

The singularities discovered here are generalised sud¬ 
den singularities: the divergence in the scale factor occurs 
at no lower order than the third derivative, a, (as in Ref. 
[9]). Note that while we can hnd potentials leading to 
these very weak singularities, and hnd potentials in the 
same family which admit big crunch singularities, it is 
not so easy to construct scalar-held models which admit 
sudden singularities with divergence in a. It is not diffi¬ 
cult to see why. At such a singularity we would require 
H to be hnite, while H diverges. This implies p oo 
on approach to the singularity, and so V(p) must diverge 
in such a way as to cancel precisely the divergence of 
p. Such potentials will likely be unphysical, and if they 
do admit sudden singularities then the sensitive cancel¬ 
lation required will likely render them a set of measure 
zero. Sudden hnite-time singularities with a —oo are in 
some sense a borderline case between singularities which 
are dehnitely inextendible (big crunch or rip) and those 
which are not, so we could view our results as indicating 
that they may not exist in general for ’reasonable’ forms 
of matter, but their generalised forms with dTa/dV^ —> oo 
for n > 3 are common. 

Although we restricted to spatially hat universes our 
arguments are unahected if the universe is open. We 
also expect hnite-time singularities to form in a closed 
universe provided they form before the universe reaches 
the expansion maximum, as in Ref. [6]. However, a mas- 
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sive scalar field (n = 2) in a closed FLRW universe can 
have quite intricate fractal behaviour that allows a subset 
of initial data to avoid a strong curvature collapse singu¬ 
larity [40], so we might expect analogously complicated 
behaviour for non-integer potentials. 

Note that we expect these singularities to be classically 
stable. This is because the formation of these singulari¬ 
ties is stable within the family of FRLW spacetimes, so 
they can only be destabilised by the formation of large 
anisotropies or inhomogeneities as ^ ^ 0. One can show 
they also generically form in the anisotropic case, and 
since sudden singularities are known to be classically sta¬ 
ble [13] we expect these even weaker singularities to be 
likewise stable. Indeed, as sudden singularities have been 
shown to be stable [14] in the sense that the homogeneous 
solution is the leading order solution of the field equa¬ 
tions containing nine independently arbitrary functions 
in the neighbourhood of the singularity we would expect 
these finite-time scalar-field singularities to be stable in 
this sense (however, the general scalar-field solution only 
contains six arbitrary spatial functions on a time slice). 
They are also probably even stable to quantum correc¬ 
tions, since the curvature on approach to the singularity 
is finite and we can construct singularities in which arbi¬ 
trarily high derivatives of the curvature are finite [19, 20]. 
This is quite unlike most other examples of exotic singu¬ 
larities known. 

A pressing question raised by this paper is what is the 
future evolution of these spacetimes beyond the point (j) = 
0? Since all polynomial curvature invariants are finite, 
the spacetime is extendible past this point. Notice that as 
^ 0 (with ^ < 0) then the scale factor approaches that 

of a universe dominated by a stiff fluid: a{t) oc t's. The 
precise form of on approach to the singular point is in 
general quite complicated and depends upon the precise 
form of V{(/)). The problem with evolving the system 


beyond (j) = 0 is that in some cases (for instance n = ^) 
the matter model breaks down beyond the point ^ = 0, 
since the naive evolution would push (j) to become strictly 
negative, which in general would lead to the expansion 
rate becoming complex. Note though that this is only the 
case for some choices of n, and there are many choices 
(e.g. n = ^) for which V{(j)) is always real. In this 
case numerical evidence suggests that iiV{(j)) is negative- 
definite for (() < 0 the spacetime collapses to a big crunch 
singularity, while if V{(j)) is positive-definite for ^ < 0 no 
such collapse occurs. Instead the universe goes through 
^ = 0 an infinite number of times. We hope to return to 
this issue in future work. Some ideas for how to approach 
this problem using distributional quantities are discussed 
in Refs. [41, 42]. 

A further question is how this model of inflation could 
be distinguished from a similar large-held model? Also, 
would these singularities reported here survive the de¬ 
tailed physics of reheating, which usually involves cou¬ 
pling the inhaton to one or more particles? This is of 
course impossible to answer without a specihc model of 
reheating, but there is no reason to think reheating would 
alter them. One way to see this is that reheating can usu¬ 
ally be modeled by adding a friction-like term, F((), into 
Eq. (3), and it is not difficult to see this should not 
generally affect the ability of the held to reach cj) = 0 
[43]. The prospects of detecting an observable signature 
of this behaviour face the same challenges as connecting 
the details of reheating to CMB observations. 
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